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HIGHLIGHTS 


•  The  MA-ES  provides  comprehensive  valuations  on  ORC  used  for  waste  heat  utilization. 

•  The  MA-ES  covers  energetic,  exergetic  and  economic  evaluations  of  typical  ORCs. 

•  The  MA-ES  is  a  general  assessing  method  without  restriction  to  specific  ORC  condition. 

•  Two  ORC  cases  of  ICE  waste-heat-recovery  are  exemplified  applying  the  MA-ES. 
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A  Multi-Approach  Evaluation  System  (MA-ES)  is  established  in  this  paper  providing  comprehensive 
evaluations  on  Organic  Rankine  Cycles  (ORC)  used  for  waste  heat  utilization.  The  MA-ES  covers  three 
main  aspects  of  typical  ORC  performance:  basic  evaluations  of  energy  distribution  and  system  efficiency 
based  on  the  1st  law  of  thermodynamics;  evaluations  of  exergy  distribution  and  exergy  efficiency  based 
on  the  2nd  law  of  thermodynamics;  economic  evaluations  based  on  calculations  of  equipment  capacity, 
investment  and  cost  recovery.  The  MA-ES  is  reasonably  organized  aiming  at  providing  a  general  method 
of  ORC  performance  assessment,  without  restrictions  to  system  configurations,  operation  modes, 
applications,  working  fluid  types,  equipment  conditions,  process  parameters  and  so  on.  Two  ORC  cases 
of  internal  combustion  engines’  (ICEs)  waste-heat-recovery  are  exemplified  to  illustrate  the  applications 
of  the  evaluation  system.  The  results  clearly  revealed  the  performance  comparisons  among  ORC  config¬ 
urations  and  working  fluids  referred.  The  comparisons  will  provide  credible  guidance  for  ORC  design, 
equipment  selection  and  system  construction. 

©  2014  Elsevier  Ltd.  All  rights  reserved. 


1.  Introduction 

The  issue  of  global  energy  inadequacy,  together  with  climate 
and  economic  issues,  has  become  more  and  more  urgent  currently. 
Countries  around  the  world  are  developing  policies  and  techniques 
in  industries  for  better  utilization  of  residual  energy  that  accounts 
a  great  part  of  original  fuel  energy.  Residual  energy  in  industries 
mainly  refer  to  low  and  medium  grade  waste  heat,  which  accounts 
for  more  than  half  of  the  total  heat  generated  [1],  and  it  widely 
exists  in  almost  every  industrial  field,  especially  in  cement  indus¬ 
try,  oil  industry,  engines,  gas  turbines,  solar  field,  biomass  and 
geothermal  fields  [2,3]. 
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Among  various  techniques  of  utilizing  residual  energy,  Organic 
Rankine  Cycle  (ORC)  is  receiving  more  and  more  attention  for  its 
high  efficiency  and  flexibility,  and  good  matching  between  work¬ 
ing  fluids  and  low  and  medium  grade  heat  sources.  Comparing 
with  steam  rankine  cycle  (RC)  under  the  same  residual  conditions, 
ORC  has  the  advantages  of  higher  thermal  efficiency,  more  net 
power,  smaller  system  volume  and  weight  as  well  as  less  complex 
expander  devices,  etc.  [4],  It  can  be  concluded  from  literature  sur¬ 
vey  that  plenty  of  theoretical  researches  on  ORC  systems  used  in 
low  and  medium  grade  heat  utilization  have  been  conducted 
focusing  on  various  applications,  working  fluid  types  [5,6],  opera¬ 
tion  modes  [27],  system  configurations  [28],  specific  components 
[7]  and  performance  parameters.  Actually,  a  reasonable  and  com¬ 
prehensive  evaluation  system  is  essential  to  ORC  systems  used  in 
waste  heat  utilization  to  provide  guiding  assessments  and  compar¬ 
isons  of  different  ORC  operations.  Among  current  evaluation 
researches,  three  main  analysis  aspects  of  ORC  performance  are 
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Nomenclature 

Abbreviations 

9 

the  heat  flux  [W/m2] 

ORC 

Organic  Rankine  Cycle 

M 

molar  mass  of  working  fluid  [kg/kmol]. 

MA-ES 

Multi-Approach  Evaluation  System 

K 

tHTC;  the  total  heat  transfer  coefficient  of  a  heat 

BSFC 

brake  specific  fuel  consumption 

exchanger  [kW/(m2  K)] 

CHP 

combined  heating  and  power  system 

A 

heat  transfer  area  [m2] 

EIP 

expander-inlet-pressure  [MPa] 

As 

heat  transfer  area  per  net  power  [m2/kW] 

PPTD 

Pinch  Point  Temperature  Difference  method 

A  Tm 

logarithmic  mean  temperature  difference 

HTC 

heat  transfer  coefficient  [W/(m2  K)] 

A„k 

annuity  of  the  investment  [Dollar] 

tHTC 

Total  heat  transfer  coefficient  [kW/(m2  K)] 

fr 

the  maintenance  and  insurance  cost  factor 

CEPCI 

Chemical  Engineering  Plant  Cost  Index 

CRF 

capital  recovery  factor 

Superscripts  and  subscripts 

EPC 

the  electricity  production  cost  [Dollar/kWh] 

0 

dead  state 

DPP 

depreciated  payback  period 

system 

the  whole  ORC  system 

NPV 

net  present  value 

1st 

the  first  law  of  thermodynamics 

SL 

the  single-loop  ORC 

2nd 

the  second  law  of  thermodynamics 

DL 

the  dual-loop  ORC 

hs 

the  heat  source 

HL 

the  high-temperature  loop  of  DL  ORC 

cs 

the  cold  source 

LL 

the  low-temperature  loop  of  DL  ORC 

f 

the  working  fluid 

hsf 

the  process  between  the  heat  source  and  working 

Symbols 

fluid 

P 

power  [kW] 

gf 

the  process  between  the  exhaust  gas  and  working 

Q 

heat  flow  [kW] 

fluid 

E 

exergy  flow  [kW] 

few 

the  process  between  the  working  fluid  and  cooling 

I 

exergy  loss  [kW] 

water 

T 

temperature  [K] 

pp 

the  pinch  point 

m 

mass  flow  rate  [kg/s] 

g 

the  exhaust  gas 

cp 

specific  heat  capacity  [kj/(kg  I<)] 

jw 

the  jacket  water 

t 

plate  thickness  [m] 

cwlc/condens  the  cooling  water  for  condensing;  the  condensing 

X 

coefficient  of  thermal  conductivity  [kW/(m  K)] 

process 

n 

efficiency  [dimensionless] 

V 

the  vapor  state 

a 

heat  transfer  coefficient  [W/(m2  K)] 

/ 

the  liquid  state 

P 

chevron  angle  [degree] 

wall/w 

the  wall-temperature  state 

P 

density  of  working  fluid  [kg/m3] 

e/evp 

the  evaporation  process 

P 

viscosity  [pa  s] 

bm 

the  bare  module 

V 

velocity  [m/s] 

exp 

the  expander 

h 

specific  enthalpy  [kj/kg] 

fp/pump 

the  working  fluid  pump 

s 

specific  entropy  [kj/(kg  K)] 

pre 

the  preheated  region  of  working  fluid 

fg 

the  friction  factor  in  exhaust  gas  side 

suph 

the  superheated  region  of  working  fluid 

Re 

Reynolds  number  [dimensionless] 

supc 

the  superheated  region  of  working  fluid  in  the  con¬ 

Pr 

Prandtl  number  [dimensionless] 

densing  process  (5-6) 

Nu 

Nusselt  number  [dimensionless] 

cond 

the  condenser/condensation 

C 

working  fluid  mass  flux  [kg/m2  s] 

cout 

the  cooling  water  flowed  out  from  the  condenser 

S 

acceleration  due  to  gravity  [m/s2] 

m 

the  process  between  the  two  working  fluids  of  two 

ifg 

the  enthalpy  of  vaporization  [kj/kg]. 

loops  in  DL  ORC 

Fr 

Froude  number  [dimensionless] 

jw/2 

the  process  between  the  jacket  water  and  working 

Co 

convection  number] dimensionless] 

fluid  in  the  low-temperature  loop  of  the  DL  ORC 

Bo 

boiling  number  [dimensionless] 

/l,expl,»7i,gl,gfl,pumpl  parameters  or  processes  of  the  high- 

X 

the  vapor  quality 

temperature  loop  of  the  DL  ORC 

Deq 

the  equivalent  diameter  of  working  fluid  channel 

f2,exp2,ri2,gX,gf2,pump2,supc2,cond2  parameters  or  processes  of 

Per 

the  critical  pressure  of  working  fluid  [MPa] 

the  low-temperature  loop 

Pe/Pmax 

the  evaporating  pressure  of  working  fluid;  EIP;  the 

of  the  DL  ORC 

maximum  pressure  of  the  ORC  cycle[MPa] 

l-13,2g,3g,5s,7s,13s  state  points 

Rp 

mean  asperity  height  [pm] 

generally  considered,  i.e.,  basic  energy  aspect,  exergy  aspect  and 
economical  aspect. 

As  to  the  basic  energy  analyses  of  ORC  systems,  Iacopo  and  Agos- 
tino  [  8  ]  described  a  specific  thermodynamic  analysis  in  order  to  effi¬ 
ciently  match  ORC  system  with  stationary  internal  combustion 
engine  (ICE).  Main  target  parameters  concerned  consisted  of  the 
ORC  efficiency,  net  power  output,  mass  flow  rate  of  working  fluid, 
and  turbine  outlet/inlet  volume  flow  ratio,  etc.  The  investigations 
were  mainly  based  on  the  basic  energy  analyses  of  ORC  systems  by 


manipulating  turbine  inlet  pressure.  Katsanos  et  al.  [9]  conducted 
a  theoretical  study  to  investigate  the  potential  improvement  of  the 
overall  efficiency  of  a  heavy-duty  truck  diesel  engine  equipped  with 
an  ORC  bottoming  system  for  recovering  heat  from  exhaust  gas. 
Based  on  basic  energy  level  analysis,  main  target  parameters  con¬ 
cerned  are  brake  specific  fuel  consumption  (BSFC)  improvement  of 
diesel,  ORC  generated  power,  working  fluid  mass  flow  rate,  cycle 
high  pressure  and  cycle  efficiency.  Zhang  et  al.  [  1 0  ]  proposed  a  novel 
system  combining  a  vehicular  light-duty  diesel  engine  with  a  dual 
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loop  ORC,  which  recovers  waste  heat  from  the  engine  exhaust, 
intake  air  and  coolant.  Main  target  parameters  include  the  overall 
net  power,  thermal  efficiency,  improvement  of  engine  effective 
power  and  improvement  of  BSFC  by  manipulating  the  engine  condi¬ 
tions,  namely,  the  heat  sources’  conditions.  Lisa  et  al.  [1 1  ]  analyzed 
and  compared  various  possible  modifications  to  simple  ORC  plant 
layout,  such  as  recuperated  cycles,  superheated  cycles,  supercritical 
cycles,  regenerative  cycles  and  their  combinations.  Based  on  basic 
energy  assessments,  six  main  thermodynamic  indexes  were  investi¬ 
gated,  including  cycle  efficiency,  specific  work,  recovery  efficiency, 
turbine  volumetric  expansion  ratio,  ORC  fluid-to-hot  source  mass 
flow  ratio  and  heat  exchangers  size  parameter.  The  investigations 
were  aimed  at  providing  useful  guidelines  for  selecting  ORC  config¬ 
uration,  appropriate  fluids  and  for  defining  operating  parameters 
with  specific  given  application. 

As  to  the  exergy  analyses  of  ORC  systems,  Huijuan  et  al.  [12] 
researched  on  transcritical  Rankine  cycles  using  refrigerant  R32 
(CH2F2)  and  carbon  dioxide  (C02)  as  the  working  fluids  for  the  con¬ 
version  of  low-grade  heat  into  mechanical  power.  Aside  from  basic 
energy  calculations,  analysis  of  the  exergy  destruction  and  loss  as 
well  as  the  exergy  efficiency  optimization  of  transcritical  Rankine 
cycle  was  conducted.  Yiping  et  al.  [13]  optimized  the  thermody¬ 
namic  performance  of  ORCs  for  waste  heat  recovery  with  exergy 
efficiency  as  an  objective  function  by  means  of  genetic  algorithm. 
One  result  showed  that  the  ORC  system  using  R236ea  has  higher 
exergy  efficiency  than  other  working  fluids  did  under  the  same 
condition  of  given  waste  heat.  Roy  et  al.  [14]  developed  a  computer 
program  to  parametrically  optimize  and  compare  the  systems  and 
second  law  efficiency,  irreversibility  of  the  system,  availability 
ratio,  work  output  and  mass  flow  rate  with  increasing  turbine  inlet 
temperature  of  ORC  under  different  conditions  of  heat  source.  They 
confirmed  that  a  non-regenerative  ORC  during  superheating  using 
R123  appears  to  be  a  good  choice  for  converting  low-grade  heat  to 
power  with  high  efficiencies,  turbine  work  output  and  low 
irreversibility.  Tchanche  et  al.  [15]  studied  on  three  modified 
engines  derived  from  simple  Rankine  engine  using  regeneration 
(incorporation  of  regenerator  or  feed-liquid  heaters)  with  exergy- 
topological  method.  Exergy  destruction  throughout  the  systems 
operating  with  R134a  was  quantified  and  illustrated  using  exergy 
diagrams.  The  results  revealed  that  greater  exergy  destruction 
occurred  in  devices  such  as  turbine,  evaporator  and  feed-liquid 
heaters  than  in  other  devices.  Similar  exergy  analyses  of  ORC 
systems  were  conducted  in  Refs.  [16-18], 

As  to  the  economical  analyses  of  ORC  systems,  Hua  et  al.  [4] 
proposed  an  ORC  model  used  to  recover  waste  heat  in  ICE  exhaust 
gas.  The  ORC  system  was  techno-economically  analyzed.  The 
screening  criteria  included  the  thermal  efficiency  of  ORC,  expan¬ 
sion  ratio,  net  power  output  per  unit  mass  flow  rate  of  hot  exhaust, 
the  ratio  of  total  heat  transfer  area  to  net  power  output  and  elec¬ 
tricity  production  cost.  Lecompte  et  al.  [19]  developed  a  thermo- 
economic  calculation  procedure  for  ORC  taking  into  account  of  part 
load  operations.  The  methodology  was  applied  to  a  CLIP  system 
driven  by  an  ICE  where  waste  heat  is  available  from  the  closed  loop 
engine  cooling  circuit.  The  economically  optimal  ORC  system  was 
firstly  sized  for  the  heat  source  and  then  was  used  in  an  hour- 
by-hour  quasi-steady  state  simulation  over  a  reference  year  to 
investigate  possible  net  energy  productions.  Schuster  et  al.  [20] 
developed  energetic  and  economic  investigations  of  ORC  systems 
used  in  various  applications.  It  revealed  that  the  efficiency  of 
waste-heat-recovery  utility  is  strongly  dependent  on  the  interest 
rate  and  the  full  load  hours  of  the  system  and  that  reliable  and  long 
operation  is  necessary  for  an  economic  operation  of  a  waste-heat- 
recovery  unit.  The  thermodynamic  and  economic  performance  of 
both  subcritical  and  transcritical  ORC  systems  were  examined  for 
low-temperature  geothermal  power  plant  by  Shengjun  et  al.  [21], 
They  conducted  a  simulation  procedure  using  five  indicators: 


thermal  efficiency,  exergy  efficiency,  recovery  efficiency,  heat 
exchanger  area  per  unit  power  output  and  levelized  energy  cost. 
The  saved  petroleum  and  reduced  C02  emission  per  year  were  also 
simply  estimated.  Gewald  et  al.  [22]  explored  rankine  cycles  as 
bottoming  waste-heat-recovery  cycles  to  enhance  electrical  effi¬ 
ciency  of  internal  combustion  engines.  Energetic,  exergetic  and 
economic  analyses  were  conducted.  The  annuities  of  capital- 
related  and  operation-related  costs  as  well  as  the  specific  costs  of 
electricity  production  were  calculated.  Pedro  and  Rogelio  [23]  have 
presented  a  simple  methodology  to  determine  the  economic,  ener¬ 
getic  and  environmental  potential  benefits  that  can  be  obtained 
from  the  implementation  of  an  MT-ORC.  Economic  indexes  like 
operational  cost,  primary  energy  consumption  savings  and 
payback  period  have  been  investigated. 

Previous  ORC  researches  as  reviewed  above  are  mostly  not 
comprehensive  enough  to  cover  all  important  aspects  of  system 
performance.  Although  each  of  the  three  ways  of  evaluations  has 
been  proposed  and  applied  to  varying  degrees  by  other  researchers 
in  many  chemical  engineering  fields,  it  is  of  great  importance  to 
select,  combine  and  complement  the  evaluations  in  order  to  form 
a  comprehensive  system  for  ORC  evaluation. 

In  this  work,  a  comprehensive  evaluation  system  of  ORC 
performance,  which  is  defined  as  a  Multi-Approach  Evaluation 
System  (MA-ES),  is  established  in  this  paper  by  combining  and 
complementing  the  existing  ORC  modeling  methods.  Based  on 
the  MA-ES,  theoretical  analyses  and  performance  comparisons  of 
ORC  systems  will  become  more  comprehensive  and  clear,  since  it 
provides  a  general  method  of  ORC  performance  assessment, 
without  restrictions  to  system  configurations,  operation  modes, 
applications,  working  fluid  types,  equipment  conditions  and  process 
parameters,  etc.  The  comparisons  will  provide  valuable  guidance 
for  ORC  design,  equipment  selection  and  system  construction. 
Finally,  two  typical  ORC  cases  of  engine  waste-heat-recovery  are 
investigated  in  this  work  to  illustrate  the  applications  of  the 
evaluation  system. 


2.  MA-ES  descriptions 

The  proposed  MA-ES  comprises  three  main  parts  covering  the 
aspects  of  ORC  performance.  They  are  illustrated  in  the  schematic 
diagram  in  Fig.  1  and  are  explained  in  detail  as  follows  one  by  one. 

Wayl  of  the  MA-ES  focuses  on  evaluations  based  on  the  1  st  law 
of  thermodynamics,  i.e.,  energy  conservation  principle.  The  aim  of 
this  part  is  to  investigate  the  energy  distributions  and  conversions 
in  quantity  level.  It  is  the  fundamental  part  for  ORC  researches  and 
is  essential  for  further  investigations.  The  most  typical  perfor¬ 
mance  indexes  in  this  part  include  net  power  Pnet  and  system 
efficiency  rjut  (the  1st  law  efficiency),  etc.  The  detailed  modeling 
is  presented  in  Section  3.2. 

Way2  of  the  MA-ES  focuses  on  evaluations  based  on  the  2nd 
law  of  thermodynamics,  which  is  to  probe  into  the  exergy  distribu¬ 
tions  of  ORC  systems.  Basic  computations  based  on  heat  balance  in 
Wayl  cannot  reflect  explicitly  the  extent  of  energy  utilization. 
Evaluations  of  system  exergy  are  indispensible  since  on  the  one 
hand,  they  can  figure  out  how  efficient  the  energy  is  converted 
and  utilized  in  quality  level  and  on  the  other  hand,  they  can  reveal 
the  exergy  losses  occurred  in  each  component  and  find  out  the 
component  with  the  maximum  exergy  loss,  pointing  out  optimiz¬ 
ing  directions.  The  typical  performance  indexes  in  this  part  include 
total  exergy  loss  fystem  and  exergy  efficiency  rj2nd  (the  2nd  law  effi¬ 
ciency),  etc.  The  detailed  modeling  is  presented  in  Section  3.3. 

Way  3  of  the  MA-ES  focuses  on  techno-economic  evaluations  of 
ORC  systems  based  on  calculations  of  components’  capacity, 
investment  and  cost  recovery.  This  part  is  significant  for  compre¬ 
hensive  evaluations  on  ORC  system  and  is  meaningful  for  system 
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Level  1  Level  2  Level  3 


Wayl 

Is1  Law  of 
Thermodynamics 
based  Evaluation 

1 _ 1 

ORCegl 

1  1 

ORC  eg2 

ORC  eg3 

1  1 

P„,  or 

Way2 

2*'  Law  of 
Thermodynamics 
based  Evaluation 

1  1 

ORCegl 

1  1 

ORC  eg2 

ORC  eg3 

*\2nd 

etc . 

Way3 

Techno-Economic 

Evaluation 

ORCegl 

1  1 

ORC  eg2 

ORCegi 

Index: 

4S. 

Cost2012 

EPC 

_ elSL. 

Fig.  1.  Schematic  diagram  of  the  ORC  Multi-Approach  Evaluation  System  (MA-ES). 


design  and  construction  by  providing  credible  predictions  and 
guidance.  Many  costs  other  than  purchased  costs  of  equipments 
are  taken  into  consideration  in  the  economic  evaluations  of  an 
integrated  ORC  system.  The  factors  associated  with  system  invest¬ 
ments  include  direct  project  expenses,  indirect  project  expenses, 
contingent  and  contractor  fees  as  well  as  auxiliary  facility 
expenses.  In  addition,  material  factors  and  pressure  factors  are 
taken  into  account  in  the  cost  evaluations.  Associating  with  ther¬ 
modynamic  calculations  in  the  first  two  ways,  the  amount  of 
investment  and  cost  recovery  period  can  be  figured  out  by  the 
techno-economic  evaluations  in  this  part.  The  detailed  modeling 
is  presented  in  Section  3.4. 

As  shown  in  the  schematic  diagram  of  Fig.  1,  three  fictitious 
examples  of  ORC  system  (ORCegl,  ORCeg2,  ORCeg3)  are  proposed 
to  better  interpret  the  evaluation  system.  Hypothetically,  they  are 
three  different  systems  running  under  the  same  given  boundary 
condition.  Evaluation  results  are  collected  in  Fig.  1  to  compare  their 
performance.  Comparing  the  three  sets  of  result,  the  value  range  of 
each  index  can  be  easily  got.  The  value  ranges  of  all  indexes  are 
graded  into  three  levels  (like  Level  1:  a1-a2,  Level  2:  a2-a3,  Level 
3:  a3-a4)  averagely  and  increasingly  in  order  to  precisely  reveal 
the  relative  performance  of  each  system.  According  to  the  1  st  law 
of  thermodynamics  based  analysis  in  the  Wayl  part,  the  ORCegl 
system  performs  poorly  with  low  Pnet  or  t/lst  which  are  ranked  in 
Levell  while  the  ORCeg3  performs  the  opposite.  According  to  the 
2nd  law  of  thermodynamics  based  analysis  in  the  Way2  part,  the 
pystem  ORCegl  is  the  lowest  which  means  the  least  exergy  loss. 
However,  if  the  t]2nd  is  chosen  as  the  representing  index  in  this  part, 
the  ORCegl  shows  the  biggest  irreversibility  with  the  lowest  rj2nd. 

Meanwhile,  the  low  level  (Level  1)  in  Way  3  evaluation  is  the 
low-cost  and  good-economy  region.  As  shown  in  this  example, 
economical  indexes  of  ORCegl  obtained  from  Way  3  evaluations 
are  all  lower  than  those  of  the  other  two  systems,  meaning  better 
economic  performance.  The  ORCeg3  system  indicates  opposite 
characteristics  with  the  largest  economic  indexes.  Overall,  the 
ORCeg2  is  a  medium  system  with  moderate  performance. 

The  evaluation  system  operates  in  this  way  for  clearly  revealing 
of  ORC  system  performance.  It  is  modeled  and  expounded  in  detail 
in  the  paragraphs  below. 

3.  MA-ES  modeling  and  analyzing 

The  MA-ES  modeling  process  in  this  part  is  introduced  based  on 
a  conventional  ORC  system.  A  basic  ORC  system  could  derive  kinds 


of  ORC  modifications,  such  as  preheated  ORCs  [8,24],  regenerative 
ORCs  [11,25]  and  dual-loop  ORCs  [10,25,26],  By  making  corre¬ 
sponding  adjustments  to  specific  correlations  and  coefficients, 
the  proposed  MA-ES  is  applicable  to  all  ORC  modified  configura¬ 
tions  because  they  originate  from  the  same  ORC  principle.  It  pro¬ 
vides  a  general  method  of  ORC  performance  assessment,  without 
restrictions  to  system  configurations,  operation  modes,  applica¬ 
tions,  working  fluid  types,  equipment  conditions  and  operation 
parameters,  etc. 

3.1.  ORC  descriptions 

The  basic  ORC  principle  is  illustrated  in  Figs.  2-4.  In  a  word,  a 
complete  ORC  system  basically  comprises  a  compressing  process 
in  the  pump  for  pressure  lifting,  a  heating  process  by  hot  sources 
for  temperature  increasing,  an  expanding  process  in  the  expander 
for  power  output  and  a  condensing  process  in  the  condenser  for 
fluid  cooling  down.  Cycles  run  continuously  in  this  way  to  generate 
continuous  power.  ORC  operations  are  typically  divided  into  two 
operating  modes  by  expander-inlet-pressure  (EIP).  They  are  sub- 
critical  systems  with  EIP  below  critical  pressure  and  transcritical 
systems  with  EIP  above  critical  pressure  [27],  as  shown  in  Figs.  3 
and  4  separately. 


Fig.  2.  Schematic  diagram  of  a  typical  ORC  system. 
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Before  establishing  the  MA-ES  models,  assumptions  are  made  as 
following:  each  component  is  considered  as  a  steady-state  and 
steady-flow  system;  the  kinetic  and  potential  energies  as  well  as 
heat  and  friction  losses  are  neglected;  the  isentropic  efficiencies 
of  the  pump  and  the  expander  are  set  at  constant  values.  Besides, 
boundary  conditions  for  simulations  should  be  set  clearly  accord¬ 
ing  to  certain  ORC  applications,  including  environmental  condi¬ 
tions,  isentropic  efficiencies  of  the  pump  and  the  expander,  heat 
and  cold  source  conditions,  crucial  operating  parameters,  pinch 
temperatures  of  heat  exchangers,  and  so  on. 


3.2.  MA-ES  Way  J;  The  1st  law  of  thermodynamics  analysis 


Wayl  shows  the  first  way  of  performance  evaluating  based  on 
energy  level  analysis.  In  this  part,  the  energy  balance  equations 
of  system  components  and  the  system  efficiency  equation  are 
listed  in  the  second  column  in  Table  1.  It  can  be  revealed  that, 
the  energy  conversion  in  each  component  is  an  arithmetic  product 
of  the  mass  flow  rate  of  worldng  fluid  and  the  enthalpy  difference 
of  inlet  and  outlet  fluid.  So  the  fundamental  calculation  target  in 
this  part  is  the  mass  flow  rate  of  the  working  fluid,  which  deter¬ 
mines  the  thermodynamic  state  of  every  point  in  the  ORC  system 
and  the  performance  of  the  whole  system.  A  common  method  for 
mass  flow  rate  calculation  and  thermodynamic  states’  computa¬ 
tion  is  the  Pinch  Point  Temperature  Difference  (PPTD)  method 
[4,8],  The  PPTD  method  was  illustrated  for  subcritical  ORC  in  pre¬ 
vious  work  by  Tian  of  our  group  [4],  In  this  work,  the  method  is  fur¬ 
ther  introduced  for  transcritical  system  since  the  modeling  process 
is  quite  different  in  the  two  system  modes. 

Take  the  transcritical  cycle  in  Fig.  4  for  instance,  the  heat 
exchange  between  the  heat  source  and  working  fluid  can  be  eval¬ 
uated  as  following  by  applying  PPTD  method.  Parameters  like  the 
temperature  of  the  inlet  heat  source  Thsin,  mass  flow  rate  of  heat 
source  rhte,  working  fluid  temperatures  at  point  1  and  point  4  have 
been  known  prior  to  this  calculation.  The  pinch  point  temperature 
of  this  heat  exchange  process  is  expected  to  be  Tpp. 

Step  1,  the  temperature  of  the  outlet  heat  source  is  assumed  as 
Thsout  (a  constant  value),  thus  the  temperature  difference  (AThs)  of 
the  inlet  and  outlet  heat  sources  is  got. 

Step  2,  the  mass  flow  rate  of  working  fluid  ( mf0 )  can  be  prelim¬ 
inarily  determined  by  the  following  Eq.  (1). 


■  CPhsItlhs (Thsin  ~  Thsout ) 
CivfT.-T.) 


(1) 


Step  3,  the  temperature  difference  (A Ths)  of  the  heat  source  is 
evenly  divided  into  N  (N  is  a  constant  value)  segments  from  T/1S[0] 
to  ThS[n].  Namely,  Thsin  is  equal  to  Ths[o],  and  Thsout  is  equal  to  T^,,]. 
Correspondingly,  the  heating  process  of  working  fluid  from  point 
1  to  point  4  is  divided  evenly  from  T^o]  to  Tfln],  and  T4  is  equal  to 
7}[o],  Tt  is  equal  to  7}[n].  Besides,  T^,]  reflects  a  random  point 
between  the  inlet  and  outlet  heat  sources.  The  step-length  is  set 
to  be  A, 


(2) 


Step  4,  the  working  fluid  temperature  at  point  /i  (7}xi])  will  be 
acquired  easily  based  on  heat  balance  with  known  parameters 
7fcr[0],  Ths[i],  rhhs,  T4  and  mf0,  and  the  temperature  difference  between 
heat  source  and  working  fluid  in  the  first  segment  is  acquired  as 

ATju/ji]  =  -  1/(1]  (3) 


Similarly,  calculations  are  conducted  gradually  in  all  the  N  seg¬ 
ments,  and  the  temperature  differences  at  any  points  (AT^,])  can 
be  obtained. 


Thermodynamic  models  usei 

d  in  Wayl  and  Way2  evaluations  of  MA-ES. 

Energy  distributions 

Exergy  losses  distributions 

Expander 

Work  fluid  pump 
Evaporator/heater 
Condenser 

System 

Efficiency 

Plunder  =  "If  ■  (Jl4  -  h5) 

Ppump  =  rhf  ■  (h,  -  h7) 

Q.heat  =  thf  •  (h4  —  h i )  Or  Qki,  =  CPhs  ■  thhs  •  (Thsin  —  Thsout) 

Qcondens  =  %  '  (h5  -  h?)  Or  Qpondem  =  Cp^,  ■  J71™,  ■  ( fW  -  hd„ ) 

'expander  =  mf  ■  (S5  -  S4)  ■  T0 

Ipump  =  rhf  ■  (s,  -  s7)  ■  To 

' heat  =  (Tgin  +  ^l)  ("gout +  84) 

hardens  =  (Ecwin  +  85)  ~  (8cwout +  87)  And  Tfjf  soun:e  =  Epwout  —  Ef.win 

n  l»=JWQr,-100% 

hind  =  Pn *  100% 
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Step  5,  the  minimum  temperature  difference  can  now  be  known 
as  AThs/jnin  by  comparing  all  the  calculated  temperature 
differences  from  AThs:fll]  to  AT^nj  (N+  1  values).  Comparing  the 
AThS/-min  and  the  expected  value  Tpp.  If  the  deviation  is  acceptable, 
it  means  the  assumed  value  of  Thsout  in  Step  1  is  reasonable.  Other¬ 
wise,  the  value  of  Thsout  should  be  modified  to  a  new  one.  The 
calculation  is  then  needed  to  be  repeated  from  step  1  to  step  4  until 
a  satisfactory  deviation  occurs.  During  the  last  iteration,  the  mass 
flow  rate  %  and  outlet  temperature  of  heat  source  ThsolU  are  suc¬ 
cessfully  determined.  The  state  parameters  of  each  point  as  well 
as  energy  distributions  and  efficiencies  are  then  obtained  accord¬ 
ing  to  Table  1 . 

3.3.  MA-ES  Way2:  The  2nd  law  of  thermodynamics  analysis 

Way2  shows  the  second  way  to  evaluate  the  system  perfor¬ 
mance  based  on  exergy  level  analysis  [29],  In  this  part,  the  system’s 
exergy  performance  can  be  evaluated  as  following. 

As  to  the  whole  system,  the  exergy  balance  equation  and  exergy 
destruction  equation  are  defined  as, 

+  Ppump  =  Pexpander  +  P*™  (4) 

jsys‘em  =  jcold-source  +  +  /pump  +  fexpa„der  (5) 

Ehot-soma  js  ^  net  jnpUt  exergy  from  the  heat  source.  Ppump  and 
Pexpander  are  the  pump  power  and  expander  power  respectively.  Isys~ 
tem  is  the  system’s  total  exergy  loss,  containing  the  exergy  lost  away 
with  outlet  cold  source  the  exergy  losses  resulted  from 

non-isentropic  compression  and  expansion  ( Ipump  and  Iexpander), 
and  the  exergy  lost  during  all  heat  transfer  processes  ( Iheat-exchangers )■ 
The  Iheat-exchangers  contains  exergy  lost  in  the  vapor  generator  (Iheat) 
and  the  condenser  ( Icondens )■  Specific  equations  of  each  process  are 
listed  in  the  third  column  of  Table  1 . 

3.4.  MA-ES  Way3:  The  techno-economic  analysis 

Way3  shows  the  third  way  to  evaluate  the  system  performance 
focusing  on  economic  analysis.  The  work  in  this  part  is  mainly 
referred  from  the  book  [30],  which  is  widely  used  for  analysis,  syn¬ 
thesis  and  design  of  chemical  processes.  The  capacities  of  the 
expander  and  the  pump  refer  to  the  expansion  power  Pexpa nder 
and  the  pump  power  consumptions  Ppump,  which  have  been 
obtained  from  Wayl  evaluation.  The  capacity  of  a  heat  exchanger 
refers  to  its  heat  transfer  area,  which  is  the  first  main  target  in 
the  part  of  economic  evaluation. 

3.4.1.  Calculations  of  heat  exchangers’  capacities 

The  calculation  process  of  heat  transfer  area  includes  three 
steps  as  described  below. 

Stepl,  choose  proper  types  of  heat  exchangers  according  to 
applications.  As  to  different  fluid  types  under  different  thermody¬ 
namic  conditions,  heat  exchangers  should  be  carefully  selected  for 
the  best  heat  transferring  and  easiest  construction.  Whichever  type 
of  heat  exchanger  is  chosen,  the  calculation  process  is  similar  as 
introduced  below.  Plate  heat  exchangers  are  chosen  here  for 
instance,  based  on  the  products  bought  in  our  research  group. 
The  plate  heat  exchangers  are  popular  in  modern  power  and  pro¬ 
cess  industries  due  to  its  advantages  and  outstanding  capability 
to  recover  heat  with  extremely  small  temperature  difference,  etc. 
[31,32], 

Step  2,  choose  proper  heat  transfer  correlations.  Various  kinds 
of  heat  transfer  correlations  have  been  studied  and  published 
according  to  different  applications.  None  of  them  is  identified  as 
the  best  one  for  all  applications.  Proper  selections  of  heat  transfer 
correlations  are  essential  to  creditable  evaluation  about  the  aimed 


capacity.  The  correlations  below  are  selected  for  the  basic  ORC 
example  mentioned  in  Figs.  2-4.  Besides,  the  hot  source  is  sup¬ 
posed  to  be  exhaust  gas  from  industry,  and  the  cold  source  is 
domestic  water. 

The  hot  exhaust  gas  is  supposed  to  stay  in  single  phase,  and  its 
heat  transfer  coefficient  (HTC)  ag  is  obtained  by  the  following 
equations  [33].  The  correlation  proposed  by  Gnielinski  takes  addi¬ 
tionally  the  effect  of  wall  roughness  into  account  by  means  of  the 
friction  factor  fg. 

fg  =  (1.82  •  lgReg  -  1.5)-2  (6) 

When  Reg  >  104  and  Reg  <  5  x  106, 

Nug  =  (4/8)  •  (Reg  -  1000)  •  Prg/[12.7  •  (fg/ 8)05  •  (Pif  3  -  1)  +  1 .07] 

(7) 

When  Reg  <  104, 

N%  =  (4/8)  ■  Reg  •  Prg/[12.7  •  (4/8)05  ■  (Pr^/3  -  1)  + 1.07]  (8) 

ag  =  Nug  ■  kgIDeqg  (9) 

Wherein,  fg  is  the  friction  factor,  and  the  subscript  g  means  the  hot 
exhaust  gas. 

Correlations  chosen  for  subcritical  systems  (Fig.  3)  are  firstly 
introduced  here.  For  the  working  fluid  in  single  phase  state,  i.e., 
the  fluid  in  the  sub-cooling  zone  (point  1  to  point  2)  and  in  the 
sup-heating  zone  (point  3  to  point  4,  point  5  to  point  6),  heat  trans¬ 
fer  coefficients  (a f)  can  be  calculated  using  the  correlation  as  fol¬ 
lows  [4,34],  It  is  a  practically  useful  semi-empirical  equation  for 
heat  transfer  in  technical  plates  based  on  the  Leveque  analogy 
and  on  experimental  evidence: 

Re  =  p  ■  v  ■  Deq/n  (10) 

Pr  =  cp  ■  n/k  (11) 

l//05  =  cos  j?/ (0.18  ■  tan /i+ 0.36  ■  sin/?+/0/cos  0)°'5  +  (1  -  cos$/(3.8  -/i)°'5 

(12) 

When  Re  <  2000,  f0  =  64/Re,  /,  =  579/Re  +  3.85 
When  Re  ^  2000 =  (1.8  ■  lg  Re  -  1.5)  2,f  =  39/Rea289 

Nu  =  0.122  ■  Pr1/3  ■  (p//V„)1/6(f  ■  Re2  •  sin  2/i)0374  (13) 

a  =  Nu  •  k/Deq  (14) 

Wherein,  Deq  in  Eq.  (10)  should  be  the  equivalent  diameter  of  work¬ 
ing  fluid  channel. 

Deq  =  2b  ftp  (15) 

b  is  corrugation  depth  and  tp  is  the  surface  enlargement  factor  of  a 
plate  heat  exchanger. 

The  cooling  water  side  HTC  (acw)  is  evaluated  by  the  same 
model  as  that  of  single-phase  working  fluid  above  [34], 

The  Cooper’s  pool  boiling  correlation  is  chosen  for  HTC  (ae)  cal¬ 
culation  of  the  working  fluid  evaporating  process  (point  2  to  point 

3). 

ae  =  1.5  *  55  ■  (Pe/Pcrf  A2-°-2lSm  ■  (-lg(Pe/Pcr))“a55  '  <f°'67  '  M^5 

(16) 

Wherein,  p^  means  the  critical  pressure  of  working  fluid  (MPa);  Rp 
means  the  mean  asperity  height  (pm),  which  is  the  mean  roughness 
between  peak  and  valley  at  surface.  According  to  the  supplied  plate 
type  heat  exchangers  by  a  European  company,  the  value  of  Rp  is  0.3 
in  this  calculation;  q  means  the  heat  flux  (W/m2),  which  is  evalu¬ 
ated  by  assuming  initial  wall  temperature,  and  then  using  the  iter¬ 
ative  method;  M  means  the  molar  mass  of  working  fluid  (kg/kmol). 
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In  condenser,  the  condensation  heat  transfer  coefficient  (ac)  of 


all  working  fluid  can  be  calculated  using  the  correlation  as  follows 
[35], 

a,  =  0.2092  *  (k,/Deq)  ■  Re?'78  •  Pr?33  ■  Wnwa„f 14  (1 7) 

Re,  =  G  •  DeJn 1 

(18) 

Pr,  =  cp,  ■  ji,/k, 

(19) 

Co  =  (p„M)-(l/x-l)as 

(20) 

Fr/  =  G2/(pf  g-Deq) 

(21) 

Bo  =  q/G  ■  i& 

(22) 

a  =  a,  •  (0.25  •  Co  045  •  Fr°2 

5  +  75  •  Bo0'75) 

(23) 

Wherein,  x  means  vapor  quality;  G  means  mass  flux  of  the  working 
fluid  (kg/m2  s);  g  means  the  acceleration  due  to  gravity  (m/s2);  i/g 
means  the  enthalpy  of  vaporization  (kj/kg). 

Correlations  chosen  for  transcritical  systems  (Fig.  4)  are  mostly 
the  same  as  those  of  the  subcritical  systems,  except  the  heating 
process  of  working  fluid  (point  1  to  point  4).  During  this  process, 
the  Jackson  correlation  below  is  used  to  calculate  the  HTC  ( ). 
The  Jackson  correlation  was  regarded  as  one  of  the  most  accurate 
correlations  since  97%  of  experimental  data  were  correlated  with 
the  accuracy  of  ±25%  [36], 

Nu  =  0.01 83  •  Re0  82  ■  Pr° 5  ■  ( pwall/p)° 3  ■  ( cp/cp)n  (24) 

n  =  0.4,T<Tw<Tpc,  1.2*Tpc<T<Tw 

n  =  0.4  +  0.2  * [(Tw/Tpc)  - 1  ] ,  T<Tpc<Tw 

n  =  0.4  +  0.2  *  [(Tw/TPc)  -  1]  *  [1  -  5  *  ( T/Tpc  -  1)],  Tpc  <  T  <  1 .2  *  Tpc,T  <  Tw 

Calculations  of  a,  Re  and  Prare  the  same  as  Eqs.  (10)  and  (11). 
Besides,  pwau  is  the  density  of  the  working  fluid  of  the  wall-temper¬ 
ature.  cp  is  the  averaged  over  cross-section  specific  heat  under 
constant  pressure,  cp=(hw-  h)/(Tw  -  T) 

Step  3,  calculate  total  heat  transfer  coefficient  (tHTC).  Based  on 
the  HTCs  of  each  side  calculated  above  and  the  heat  conductivities 
of  heat  exchanger  materials,  tHTC  ( K )  can  be  calculated  as, 

1  /K  =  1 000/ahot_side  +  t/X  + 1 000/OcoKWfe  (25) 

Step  4,  the  capacity  of  a  heat  exchanger,  namely,  the  heat  transfer 
area,  is  obtained  as 

A  =  Q_/K/ATm  (26) 

In  addition,  heat  transfer  area  per  net  power  A,  is  an  index  reflecting 
the  system’s  compactness.  It  influences  the  system’s  economic  per¬ 
formance  to  some  extent  [37], 

AS=A/Pnet  (27) 

3.4.2.  Investment  calculations 

Knowing  the  capacity  of  each  component  of  an  ORC  system, 
investment  modeling  is  conducted  in  three  steps  below  [30]: 

Step  1 ,  calculations  of  ORC  capital  cost. 

The  capital  cost  of  an  integrated  chemical  plant  takes  into  con¬ 
sideration  many  costs  other  than  the  purchased  cost  of  the  equip¬ 
ment.  Factors  associated  with  capital  cost  includes  direct  project 
expenses  (such  as  equipment  free-on-board  cost,  material  required 
for  installation  and  labor  to  install  the  system),  indirect  project 
expenses  (such  as  freight,  insurance,  taxes,  construction  overhead 
and  contractor  engineering  expenses),  contingency  and  contractor 
fee,  and  auxiliary  facilities  (such  as  site  development,  auxiliary 
building,  off-sites  and  utilities).  In  addition,  the  material  factors 


and  pressure  factors  are  taken  into  account  in  the  cost  evaluation. 
The  investment  model  of  each  component  as  well  as  coefficients 
referred  is  listed  in  Table  2  referring  to  [30],  The  coefficients  are 
selected  based  on  the  above  capacity  evaluation  of  each  compo¬ 
nent.  Besides,  the  components’  types  are  required  in  order  to 
choose  proper  coefficients  from  the  data  base  in  [30]  which  was 
obtained  from  a  survey  of  equipment  manufacturers.  Therefore, 
the  coefficients  listed  in  Table  2  are  mainly  chosen  for  the  two  case 
studies  in  Section  4  (MA-ES  Applications). 

The  capital  cost  based  on  economic  situation  in  year  2012  is 
obtained  below, 

COSt2001  =  Cjmex p  +  Cbmjp  +  Cbm,preh  +  Cfim.e  +  Cbm^uph  +  Cbm.c  (28) 

Cost20i2  =  Cost20oi  •  CEPChon/CEPChom  (29) 

Wherein,  CEPCI200i  =  397,  CEPCI2012  =  584.6  (CEPCI  means  Chemical 
Engineering  Plant  Cost  Index).Step  2,  calculation  of  annual  invest¬ 
ment  on  ORC.Capital  recovery  factor  (CRF)  is  estimated  based  on 
the  following  equation. 

CRF  =  i  ■  (1  +  j)tfme/((i  _|_  j)nme  _  i)  (30) 

The  annuity  of  the  investment  Ank  can  be  obtained  using  equation 
Ank  =  Cost2 012  •  CRF  (31) 

i  is  interest  rate,  whose  value  is  5%;  time  is  economic  life  time, 
whose  value  is  15  years  [20]; 

Step  3,  EPC  calculation.  Electricity  production  cost  (EPC)  can  be 
obtained  by  equation 

EPC  =  (Ank  +fK  ■  COSt20u)/((PeXp  -  Pjp)  ■  hfuU-load)  (32) 

fK  is  operation,  maintenance  and  insurance  cost  factor,  whose  value 
is  1.65%;  hfuii-ioad  is  full-load  operation  hours,  whose  value  is  7500  h. 

Step  4,  DPP  calculation.  The  depreciated  payback  period  (DPP)  is 
the  number  of  time  periods  (years)  that  are  required  to  recover  the 
initial  investment  with  all  cash  flows  discounted  back  to  the  start¬ 
up  year  [40,41  ].  It  is  given  by 

DPP  =  -  ln(l  -  k  ■  Cost2m2/Fn0 ) /  In  (1  +  k)  (33) 

F„ o  =  Ep  ■  ((Pexp  -Pfp)-  hjm oad)  ~  h  •  COSt20 12  (34) 

The  discount  rate  k  is  5%  in  this  text;  Cost20 12  is  the  initial  invest¬ 
ment;  Fno  is  the  net  cash  flow  for  year  n  without  considering  the 
inflation  rate  and  is  assumed  to  remain  constant  during  the  lifetime. 
In  an  ORC  system,  the  net  cash  flow  Fn0  is  the  difference  value 
between  benefit  and  cost  annually.  The  benefit  specifically  repre¬ 
sents  the  revenue  from  selling  produced  electricity  at  a  price  of 
0.2$/kW  h,  and  the  annual  cost  refers  to  the  operation,  maintenance 
and  insurance  cost  which  is  set  asfK  ■  Cost2o  12  according  to  Eq.  (32). 

Step  5,  NPV  calculation.  The  net  present  value  (NPV)  is  another 
frequently-used  economic  index  that  sums  the  discounted  cash 
flows  of  an  energy  system.  It  integrates  and  converts  at  the  same 
time  amounts  of  money  of  various  time  periods  [42],  It  is  given  by 

NPV  =  -Cost20n  +  J2Fn/(  1  +  k)n  (35) 


Fn  =  Fno  ■  (1  +  r)n  (36) 

n  is  the  time  period  (year);  N  is  the  number  of  years  for  which  the 
economic  evaluation  is  requested;  F„  is  the  net  cash  flow  for  year  n 
considering  the  inflation  rate  r  (2%).  A  project  investment  should  be 
realized  only  if  NPV  >  0,  and  the  best  project  should  be  the  one  with 
the  highest  NPV.  It  provides  a  basic  criterion  for  economic 
evaluation. 
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Investment  models  used  in  Way3  evaluatioi 


Expander 


lg  CPT=X,e  +  xfXPe  +  KsA  lg  P'f 


Evaporator/heater 


Condenser 


lg  Cp,ump  =  K-\p  +  K2.p  lg  Pp  +  K3,p(lg  Pp)2 
Fbm.pump  =  B,.p  +  B2.p  •  F„,pump  •  P,.pump 
lg  FP,pump  =  C,.p  +  C2.p  lg  PP  +  C3p(lg  Pp)2 


+  1^2,h  lg  Ah  +  X3./,(lg  A/,)2 
'  B2Jl  ■  FM_heater 


=  K°"+  KiXX™ KlAl&A:? 

'B7jc  +  B2jCFMfondPns 


K,  =  2.705 
K2  =  1.440 
K3  =  -0.177 
It,  =  3.870 
K2  =  0.316 
K3  =  0.122 


K,  m  4.666 
K2  m  -0.156 
K3  =  0.155 
K i  =  4.666 
K2  m  -0.156 
K3  =  0.155 


Cr/C2/C3 


Ci  =  -0.245 
C2  =  0.259 
C3  =  -0.014 


Bi/B2 


/  /  6.2 


B,  =  1.89  2.35  / 

B2  =  1.35 


B,  =  0.96  2.45  / 

B2  =  1.21 


B,  =  0.96  2.45  / 

B2  =  1.21 


Other  useful  economic  indexes  could  be  extended  based  on  cap¬ 
ital  investment  above  in  this  section  in  order  to  investigate  more 
detailed  economic  performance  from  different  perspectives  [30], 

3.5.  MA-ES  evaluation  procedure 

Summing  from  Sections  3.1  to  3.4,  the  MA-ES  evaluation 
procedure  is  basically  shown  in  the  flow  diagram  of  Fig.  5.  This 
procedure  clearly  reflects  the  methodology  of  the  MA-ES  with 
ORC  performance  calculations  and  optimizations.  Besides,  the 
principle  in  this  procedure  is  applicable  to  other  ORC  modifications 
by  adjusting  specific  parameters  correspondingly  using  MA-ES. 

4.  MA-ES  applications 

4.1.  Case  Study  A:  Evaluations  and  comparisons  of  ORC  working  fluids 
and  operation  modes 

In  this  case,  a  bottoming  ORC  system  is  constructed  to  recover 
waste  heat  of  exhaust  gas  released  from  a  heavy-duty  diesel 


Fig.  5.  MA-ES  modeling  flow  diagram. 


engine.  In  a  typical  diesel  engine,  more  than  one  third  energy  of 
the  engine  fuel  is  lost  through  exhaust  gas,  so  it  is  of  great  eco¬ 
nomic  and  environmental  importance  to  recover  waste  heat  from 
the  exhaust  [38,39],  The  bottoming  ORC  system  configuration  is 
the  same  as  shown  in  Fig.  2  above,  including  two  heat  exchangers 
for  vapor  generation  and  condensation  respectively,  a  pump  and  an 
expander.  The  subcritical  and  transcritical  systems  are  both  inves¬ 
tigated  with  T-S  figures  just  the  same  as  Figs.  3  and  4.  Four  work¬ 
ing  fluids  are  selected  for  their  good  performance  as  explored  in 
the  previous  work  [4],  which  are  R123,  R245fa,  R134a  and  R32. 
Their  thermodynamic  properties  can  be  obtained  from  the 
database  of  the  software  Engineering  Equation  Solver  (EES).  The 
MA-ES  built  above  is  applied  in  this  case  for  practical  use. 

4.1.1.  System  settings  of  Case  A 

In  addition  to  basic  system  assumptions  made  above,  more 
assumptions  as  well  as  boundary  conditions  are  required  to  be 
clarified  for  the  given  application.  The  heat  source  studied  in  this 
case  is  diesel  exhaust  under  the  rated  condition.  The  composition 
of  the  exhaust  gas  on  mass  basis  has  been  calculated  as: 
C02  =  15.1%,  H20  =  5.5%,  N2  =  71.6%,  02  =  7.8%.  This  composition  is 
used  to  evaluate  gas  properties.  Detailed  parameters  setting  for 
system  operation  are  listed  in  Table  3.  Besides,  structure  dimen¬ 
sions  of  heat  exchangers  that  derived  from  series  products  are 
listed  in  Table  4.  The  EIP  of  the  working  fluid  is  selected  as  the 
manipulated  variable.  In  subcritical  cycles,  it  is  increased  from 
the  pressure  slightly  higher  than  the  condensing  pressure  to  the 
pressure  closing  to  the  critical  point.  In  transcritical  cycles,  it  is 
increased  from  pressure  slightly  above  the  critical  point  to  10  MPa. 

4.1.2.  Evaluation  results  of  Case  A 

Evaluation  results  about  six  main  indexes  of  Case  Study  A  are 
summarized  and  compared  in  Figs.  6  and  7.  Moreover,  Appendix 
A  complements  important  results  of  this  case  study  for  further 
understanding  about  the  evaluation  process.  Appendix  A  com¬ 
prises  calculation  parameters  of  four  fluids  under  both  subcritical 
and  transcritical  modes,  including  mass  flow  rate  (rii/),  enthalpy 
values  of  key  points  (ht-h7),  work  capacity  of  key  equipments 
(PMp,  P/p),  energy  absorbed  from  heat  source  (qa),  system  exergy 
loss  (P-ystt'm),  exergy  loss  ratios  of  components  ( RatioIvg ,  RatioIexp, 
Ratiolfp,  Ratiolconll  and  Ratiolcout),  heat  transfer  coefficients  (K),  heat 
transfer  areas  (A)  and  the  annuity  of  the  investment  (Anp).  Param¬ 
eters  of  operations  with  EIP  of  3  MPa  and  10  MPa  are  chosen  to 
represent  subcritical  and  transcritical  performance  separately. 
Combining  Figs.  6  and  7  and  Appendix  A,  the  evaluation  results 
of  Case  Study  A  can  be  concluded  as  follows. 
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conditions  of  Case  Study  A. 


Parameters  Values 

Hot  source:  exhaust  temperature  [K]  792.15 

Hot  source:  exhaust  pressure  [MPa]  0.12 

Hot  source:  exhaust  mass  flow  [kg/s]  0.29 

Cold  source:  water  temperature  [K]  298.15 

Cold  source:  water  mass  flow  [kg/s]  4 

TDPP  [K]  30 


Parameters 

Expander  efficiency  [/] 

Pump  efficiency  [/] 

Condensing  temperature  [K] 

Expander  inlet  temperature  E1T  [K] 
Ambient  temperature  (dead  state)  [K] 
Ambient  pressure  (dead  state)  [MPa] 


500 

298.15 

0.101 


Table  4 

Main  geometric  dimensions  of  plate  heat  exchangers. 
Parameters 

Chevron  angle,  p  [degree] 

Surface  enlargement  factor,  <p  [/] 

Plate  width,  Lw  [m] 

Plate  thickness,  t  [m] 

Mean  asperity  height,  Rp  [pm] 

Corrugation  depth,  b  [m] 

Equivalent  diameter  of  liquid  side,  Deq  [m] 
Equivalent  diameter  of  gas  side,  Deqg  [m] 
Coefficient  of  thermal  conductivity,  X  [kW/(mK)[ 


60 

1.117 

0.119 

0.0003 

0.3 

0.00224 

0.004 

0.009 

0.00163 


Level  1  Level  2  Level  3 


Fig.  7.  Evaluations  of  transcritical  ORC  systems  of  Case  Study  A. 


The  results  are  graded  into  three  levels  by  comparing  the 
performance  of  four  working  fluids  so  as  to  clearly  compare  their 
relative  merits.  In  subcritical  operations  shown  in  Fig.  6,  R123 
and  R245fa  gain  clearly  better  performance  based  on  the  1st  law 
of  thermodynamics  analysis,  with  higher  net  power  output  and 
system  efficiency  ranking  across  Level  2  to  Level  3.  These  advanta¬ 
ges  are  more  obvious  in  transcritical  operations  as  shown  in  the 
first  two  rows  of  Fig.  7. 

Similarly,  the  exergy  efficiency  of  R123  is  much  higher  than 
those  of  the  remaining  fluids  under  both  subcritical  and  transcrit¬ 
ical  conditions.  It  is  ranked  the  highest  in  Level  3,  with  maximal 
value  of  38.5%  with  EIP  of  10  MPa.  The  R245fa  gets  the  second  best 
performance  of  exergy  efficiency.  Moreover,  the  exergy  is  mainly 
lost  during  the  heat  transfer  between  the  exhaust  gas  and  the 
worldng  fluid,  as  shown  in  Appendix  A,  and  similar  results  were 
proved  by  other  scholars  [29],  Under  subcritical  conditions,  the 
second  greatest  exergy  loss  occurs  in  condensation  process  during 
the  heat  transfer  between  the  working  fluid  and  the  cooling  water. 
However,  the  exergy  loss  in  the  expansion  process  exceeds  that  of 
condensation  process  when  the  EIP  increases  from  subcritical  pres¬ 
sure  to  transcritical  pressure. 

The  evaluations  of  heat  transfer  areas  per  net  power  As  are 
shown  in  the  fourth  rows  in  Figs.  6  and  7.  Smaller  As  means  less 
volume  and  weight  of  the  system  and  better  compactness.  Besides, 
since  heat  exchangers  are  key  components  of  an  ORC  system, 
smaller  heat  areas  are  beneficial  to  system  cost  if  manufacturing 
complexity  is  ignored.  As  revealed  in  Fig.  6  and  7,  As  of  four  fluids 
show  more  or  less  the  opposite  trends  to  the  first  three  indexes. 
Namely,  the  As  values  of  R1 23  and  R245fa  are  all  at  very  low  level 
under  both  subcritical  and  transcritical  modes.  The  first  reason  is 
that  they  get  more  net  power  as  shown  in  the  first  rows.  The 
second  reason  is  that  system  heat  is  transferred  more  efficiently 
with  higher  transfer  coefficients,  as  shown  in  Appendix  A.  According 
to  Eq.  (27),  the  areas  per  power  are  thus  apparently  lower. 

The  electricity  production  costs  (EPC)  of  R123  and  R245fa  are 
much  lower  than  those  of  the  other  two  working  fluids  in  both  sub¬ 
critical  and  transcritical  operation  modes.  They  are  all  ranked  at 
Level  1  that  represents  a  better  economic  level.  However,  the 
annuities  of  the  investments  (Ank)  of  R123  and  R245fa,  which  rep¬ 
resent  system  investments,  are  higher  than  those  of  the  other  two 
working  fluids,  as  indicated  in  Appendix  A.  According  to  Eq.  (32), 
the  low  electricity  production  costs  are  result  from  high  net  power. 

According  to  the  DPP  comparison  in  the  last  row  of  Fig.  6,  R123 
reveals  the  absolutely  best  economic  performance  since  it  requires 
the  shortest  payback  time  which  ranges  from  6.4  to  7.6  years. 
Besides,  the  system  needs  10.1-13.3  years  to  recover  the  initial 
investment  applying  R245fa.  It  is  actually  a  quite  long  period  for 
a  real  project  but  it  is  still  a  first-level-ranking  one  comparing  to 
that  applying  R134a  and  R32.  Additionally,  the  NPV  in  Appendix 
A  further  reveals  the  economic  validity  of  each  kind  of  fluid  in  sub¬ 
critical  mode.  According  to  the  NPV  definition,  the  project  could 
not  be  economically  realized  applying  R32  when  EIP  is  3  MPa, 
and  R123  is  proved  to  be  the  best  one  with  the  highest  NPV. 

However,  the  economic  performance  of  R134a  and  R32 
improves  great  when  the  system  turns  to  transcritical  mode. 
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Besides  of  the  obviously  decreased  EPCs,  their  DPPs  decrease 
greatly  from  15-45  years  level  to  13-15  years  level.  It  could  be 
explained  by  the  increasing  net  power  and  the  decreasing  specific 
area  which  together  reduce  the  investment  and  raise  the  benefits. 
However,  the  high  EIP  in  transcritical  mode  does  not  obviously 
benefit  the  investment  recovery  for  R123  and  R245fa.  The  reason 
is  that  for  these  two  fluids,  the  benefits  brought  by  the  increasing 
new  power  and  decreasing  specific  area  do  not  exceed  much  the 
costs  brought  by  the  bigger  capacities  of  the  expander  and  pump 
which  result  in  higher  investment.  However  looking  from  the 
NPV  perspective,  all  four  systems  could  be  economically  realized 
under  transcritical  pressure  of  10  MPa,  and  R245fa  and  R123  are 
proved  to  be  better  alternatives  wither  higher  NPV. 

In  conclusion,  four  conventional  different  working  fluids  and 
two  typical  operation  modes  are  evaluated  and  compared  in  this 
case  study  applying  the  proposed  MA-ES.  R123  and  R245fa  are 
evaluated  as  the  most  appropriate  working  fluids  for  this  applica¬ 
tion  with  excellent  performance  in  three  aspects  of  ORC  system.  In 
addition,  it  can  be  concluded  by  comparing  the  two  result  figures 
(Figs.  6  and  7)  and  data  in  Appendix  A  that,  performance  of  trans¬ 
critical  cycle  is  mostly  better  than  that  of  subcritical  cycle  for  a  cer¬ 
tain  fluid  under  the  given  heat  source  condition.  However,  it 
should  be  noted  that  if  the  subcritical  and  transcritial  modes  need 
to  be  investigated  in  practical  applications,  a  single  map  could  be 
much  clearer  with  the  same  scale  by  combining  Figs.  6  and  7. 

4.2.  Case  Study  B:  Evaluations  and  comparisons  of  ORC  configurations 

In  a  typical  diesel  engine,  less  than  45%  of  fuel  energy  might  be 
converted  into  useful  work  output  from  crankshaft,  and  the 
remaining  energy  is  mainly  lost  through  exhaust  gas  and  jacket 
water.  However,  the  energy  in  the  two  heat  sources  differs  greatly 
in  temperature  regime,  namely,  the  energy  quality.  The  tempera¬ 
ture  of  exhaust  gas  is  about  500-800  K,  while  that  of  jacket  water 
is  below  400  K.  Two  configurations  of  bottoming  ORC  system  -  Sin¬ 
gle  loop  ORC  and  Dual  loop  ORC  -  have  been  explored  for  heat 
recovery  from  both  heat  sources.  Single  loop  ORC  systems  were 
discussed  in  Refs.  [8,29]  using  jacket  water  to  preheat  working 
fluid  before  being  evaporated  by  exhaust  gas.  Dual-loop  ORC  sys¬ 
tems  were  researched  in  Refs.  [10,25,26,28],  Waste  heat  from 
two  heat  sources  were  recovered  using  cascade  loops  according 
to  their  temperature  gradients.  The  two  ORC  strategies  showed 
respective  advantages  in  engine  applications  by  assessing  only  a 
few  performance  indexes.  In  this  case,  their  comprehensive  perfor¬ 
mance  can  be  assessed  and  compared  applying  the  proposed 
evaluation  system  to  further  demonstrate  the  application  values 
of  the  MA-ES. 

In  this  case,  two  bottoming  ORC  systems,  a  single-loop  (SL)  ORC 
and  a  dual-loop  (DL)  ORC,  are  constructed  to  recover  waste  heat 
from  both  exhaust  gas  and  jacket  water  released  from  the  same 
heavy-duty  diesel  engine  as  used  in  Case  A.  Figs.  8  and  9  show 
the  schematic  diagrams  the  two  ORC  systems  proposed.  In  the  sin¬ 
gle-loop  ORC  system,  working  fluid  is  firstly  preheated  by  jacket 
water  and  then  evaporated  and  superheated  by  exhaust  gas  within 
a  single  ORC  circuit,  as  shown  in  Fig.  8.  The  dual-loop  ORC  system 
contains  a  high-temperature  loop  (HL)  ORC  and  a  low-temperature 
loop  (LL)  ORC,  as  shown  in  Fig.  9.  The  HL  ORC  is  designed  to  recover 
the  high-temperature  part  of  exhaust  waste  heat,  which  guaran¬ 
tees  the  system’s  safety.  The  LL  ORC  is  designed  to  recover  the 
engine  waste  heat  in  jacket  water,  residual  heat  from  the  HL  ORC 
during  its  condensation  process  and  low-temperature  part  of  the 
exhaust  waste  heat  successively. 

4.2.1.  System  settings  of  Case  B 

Boundary  conditions  in  this  case  study  such  as  ambient  condi¬ 
tions,  exhaust  gas  conditions  and  heat  exchangers’  structures  are 


the  same  as  in  Case  A.  Additional  settings  for  the  two  ORC  opera¬ 
tions  in  this  case  are  listed  in  Table  5.  In  the  single-loop  (SL)  ORC 
system,  R245fa  is  adopted  as  the  working  fluid  for  its  excellent  per¬ 
formance  revealed  in  Case  Study  A.  The  subcritical  and  transcritical 
operations  are  both  considered  based  on  this  ORC  configuration.  In 
subcritical  cycles,  EIP  is  increased  from  1  MPa  to  the  pressure  clos¬ 
ing  to  the  critical  point.  In  transcritical  cycles,  EIP  is  increased  from 
the  pressure  slightly  above  the  critical  point  to  1 0  MPa.  In  the  dual¬ 
loop  (DL)  ORC  system,  water  and  R245fa  are  adopted  as  the  work¬ 
ing  fluid  for  high  loop  and  low  loop  respectively  based  on  investi¬ 
gations  of  fluids  screening  [25],  Both  loops  are  supposed  to  be 
operated  under  subcritical  pressures.  EIP  of  the  high  loop  increases 
from  1  MPa  to  6  MPa,  since  the  pressures  above  6  MPa  will  result 
that  the  working  vapor  after  expansion  falls  into  two-phase  region. 
EIP  of  low  loop  is  set  at  0.58  MPa. 

4.2.2.  Evaluation  results  of  Case  B 

Evaluation  results  about  five  main  indexes  of  Case  Study  B  are 
summarized  and  compared  in  Fig.  10.  Appendix  B  complements 
key  results  of  this  case  study.  Overall,  the  performance  of  the  trans¬ 
critical  cycle  is  much  better  than  that  of  the  subcritical  one  consid¬ 
ering  only  the  Wayl  and  Way2  evaluations  when  single-loop  ORC 
system  is  adopted.  These  regularities  are  quite  similar  to  those 
explained  in  Case  Study  A.  However,  the  dual-loop  system  changes 
the  performance  dramatically  as  summarized  as  follows. 

Firstly,  the  dual-loop  system  can  gain  more  than  twice  the  net 
power  of  single-loop  systems.  Total  net  power  of  two  cascade  loops 
is  within  35-38  kW  range  under  the  pressure  concerned.  This  is 
quite  meaningful  for  internal  combustion  engines  that  are  stuck 
in  the  bottleneck  period  of  power  and  efficiency  lifting.  As  for  sys¬ 
tem  efficiency  it  is  ranked  at  mid-high  level  of  the  dual-loop 
system  comparing  to  that  of  the  single-loop  system.  However, 
the  dual-loop  system  shows  no  obvious  advantage  in  the  efficiency 
performance  over  the  single  loop  systems,  especially  the  transcrit¬ 
ical  one.  Besides,  the  system  efficiency  does  not  increase  greatly  as 
the  EIP  of  high  loop  increases.  It  only  changes  within  a  small  range 
about  11.5-13%  when  the  high  loop  EIP  increases  from  1  MPa  to 
6  MPa.  After  all,  the  dual-loop  system  performs  well  enough  in 
the  Wayl  evaluation  part. 

Exergy  efficiency  is  another  great  advantage  of  the  dual-loop 
system  with  a  high  value  range  of  42.1-45.0%,  which  is  right  over 
that  of  the  single-loop  system.  Total  exergy  losses  of  the  whole 
dual-loop  ORC  system  is  roughly  twice  of  that  of  the  single-loop 
system,  as  shown  in  Appendix  B.  Higher  exergy  loss  in  the  dual¬ 
loop  system  are  mainly  resulted  from  larger  amount  of  component 
needed  in  the  cycle.  But  higher  exergy  efficiencies  reveal  that  the 
waste  energy  is  more  efficiently  utilized  in  the  dual-loop  ORC 
system. 

The  As  value  range  of  dual-loop  system  is  similar  to  that  of  the 
transcritical  single-loop  system.  Both  of  them  are  ranked  at  the 
mid-high  levels,  which  are  higher  than  that  of  the  subcritical  sin¬ 
gle-loop  system. 

Although  the  amount  of  the  component’s  in  dual-loop  system  is 
definitely  larger  than  that  of  the  single-loop  system,  its  electricity 
production  costs  (EPC)  are  much  lower  with  a  range  of  about  0.67- 
0.73  Dollar/kWh. 

The  biggest  advantage  of  the  SL-Subcritical  system  is  the  short¬ 
est  depreciated  payback  time  (DPP)  it  needs.  According  to  the  last 
row  of  economic  index  comparison  in  Fig.  10,  only  about  2.8- 
4.5  years  are  needed  under  the  given  pressures  using  SL-Subcritical 
configuration.  Although  the  net  power  output  is  the  lowest  of  this 
system  configuration  comparing  to  those  of  the  other  two,  the 
components’  capacities  in  this  system  are  much  lower,  leading  to 
the  least  initial  investment  and  the  shortest  payback  time. 

On  the  contrary,  the  SL-Transcritical  system  requires  the  lon¬ 
gest  payback  time  which  is  within  18.7-22.2  years,  while  the  DL 
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system  shows  moderate  performance  concerning  the  DPP  compar¬ 
ison.  Generally  the  DPP  index  separates  the  three  configurations 
obviously  by  ranking  them  into  different  levels.  In  addition, 
according  to  Appendix  B,  NPV  values  of  the  SL-Transcritical  system 
are  both  negative  under  the  two  pressures  selected,  proving  that 
they  could  not  be  economically  realized  in  these  situations.  The 
SL-Subcritical  and  DL  systems  are  proven  to  be  valid  to  some 
extent  with  high  NPVs. 

Another  meaningful  index  is  the  heat  recovery  amount,  as 
shown  in  Appendix  B.  It  shows  that,  two  ORC  configurations  pro¬ 
posed  in  this  case  study  are  able  to  recover  almost  equal  amount 
of  waste  energy  from  the  exhaust  gas  source  (Qg).  But  the  dual-loop 
system  is  much  more  efficient  in  recovering  waste  energy  from  the 
jacket  water  source  (Q,w).  Its  recovery  amount  is  ten  times  more 
than  that  of  the  single-loop  system. 

In  conclusion,  three  ORC  configurations  are  evaluated  and  com¬ 
pared  in  this  case  study  applying  the  proposed  MA-ES.  The  dual¬ 
loop  ORC  system  is  evaluated  as  a  relatively  proper  configuration 
for  waste  heat  recovery  in  this  case  study.  It  performs  excellently 
during  Wayl  and  Way2  evaluations  with  high  thermodynamic 
indexes,  while  its  economic  performance  is  not  much  outstanding 


Table  5 

Boundary  conditions  of  Case  Study  B. 


Parameters  Values 

Temperature  of  inlet  jacket  water  [K]  356.5 

Pressure  of  inlet  jacket  water  [MPa]  0.12 

Mass  flow  of  inlet  jacket  water  [kg/s]  2.7 

Difference  between  EIT  and  inlet  exhaust  temperature  [K]  40 

TDPP  between  exhaust  gas  and  working  fluid  [K]  30 

TDPP  between  jacket  water  and  working  fluid  [I<]  5 

Condensing  temperature  of  single-loop  ORC  [K]  310 

Condensing  temperature  of  double-loop  ORC  [K]  348/310 


compared  to  single-loop  systems.  So  more  optimizations  are 
needed  and  investors  should  weigh  up  all  indexes  to  make  the 
decision. 

5.  Conclusions 

A  comprehensive  evaluation  system  of  ORC  performance  (MA¬ 
ES)  is  established,  which  provides  a  general  method  of  ORC  perfor¬ 
mance  assessment,  without  restrictions  to  system  configurations, 
operation  modes,  applications,  working  fluid  types,  equipment 


Fig.  9.  Schematic  diagram  of  the  dual-loop  (DL)  ORC  system  of  Case  Study 
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conditions  and  process  parameters,  etc.  The  evaluation  system  cov¬ 
ers  three  main  aspects  of  typical  ORC  performance:  basic  thermo¬ 
dynamic  evaluations  based  on  the  1st  law  of  thermodynamics  and 
the  2nd  law  of  thermodynamics,  economic  evaluations  based  on 
calculations  of  equipment  capacity,  investment  and  cost  recovery. 
The  evaluations  will  provide  valuable  guidance  for  ORC  design, 
equipment  selection  and  system  construction. 

Two  typical  ORC  cases  of  engine  waste-heat-recovery  are 
investigated  in  this  work  to  illustrate  the  applications  of  the 
MA-ES.  Under  the  first  case  condition,  the  performance  of  the 
transcritical  cycle  is  much  better  than  that  of  subcritical  cycle 
for  a  certain  fluid.  And  comparing  with  other  two  common 


fluids,  R123  and  R245fa  are  evaluated  as  more  appropriate  work¬ 
ing  fluids  for  this  application  with  excellent  performance  in 
three  aspects  of  ORC  system.  Under  the  second  case  condition, 
the  dual-loop  ORC  system  is  a  relatively  proper  configuration 
for  waste  heat  recovery,  which  performs  excellently  during 
thermodynamic  evaluating  processes  and  shows  moderate  per¬ 
formance  in  economic  evaluations. 
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Appendix  A 

Evaluation  results  of  Case  Study  A. 


Units 

Subcritical  ORC  (pmax 

=  3  MPa) 

Supercritical  ORC  (pma 

5=  10  MPa) 

R32 

R134a 

R245fa 

R123 

R32 

R134a 

R245fa 

R123 

rhf 

kg/s 

0.302 

0.397 

0.391 

0.476 

0.344 

0.468 

0.506 

0.669 

hr 

kj/kg 

266.1 

102.9 

248.3 

239 

275.5 

110.3 

255 

245 

h2 

kj/kg 

293 

186.5 

417.4 

404.8 

/ 

/ 

/ 

/ 

h3 

kj/kg 

508.5 

279.1 

488.8 

468.7 

/ 

/ 

/ 

/ 

h4 

kj/kg 

743.1 

465.1 

615.9 

541.5 

688.6 

414.6 

537.2 

458.4 

h5 

kj/kg 

726.5 

433.9 

566.7 

492.5 

626.8 

365.4 

487.3 

418.2 

h6 

kj/kg 

514.5 

269 

430 

404 

514.5 

269 

430 

404 

h7 

kj/kg 

265 

100.7 

245.6 

236.5 

265 

100.7 

245.7 

236.5 

Pex  p 

kW 

4.99 

12.36 

19.25 

23.27 

21.26 

23 

25.27 

26.89 

PfP 

kW 

0.336 

0.896 

1.038 

1.186 

3.607 

4.49 

4.684 

5.712 

kW 

144.1 

143.8 

143.9 

143.9 

142.1 

142.4 

142.9 

142.8 

%tem 

kW 

59.55 

52.71 

45.97 

42.1 

46.33 

45.51 

43.49 

42.89 

Ratiolgf 

% 

69.47 

61.9 

59.11 

57.79 

65.6 

60.81 

62.69 

62.37 

RatioIexp 

% 

2.235 

6.55 

12.42 

17.13 

15.28 

16.52 

20.92 

25.25 

Ratiolfp 

% 

0.109 

0.323 

0.435 

0.541 

1.474 

1.874 

2.057 

2.533 

RatioIc ond 

% 

25 

27.98 

24.68 

21.09 

14.38 

17.5 

10.97 

6.465 

Rationlcout 

% 

3.186 

3.247 

3.36 

3.446 

3.27 

3.298 

3.364 

3.374 

% 

kW/(m2  K) 

/ 

/ 

/ 

/ 

0.090 

0.089 

0.086 

0.0854 

Kgfsre 

kW/(m2  K) 

0.085 

0.085 

0.096 

0.097 

/ 

/ 

/ 

/ 

%_evp 

kW/(m2  K) 

0.107 

0.110 

0.116 

0.119 

/ 

/ 

/ 

/ 

Kgfsuph 

kW/(m2  K) 

0.109 

0.110 

0.114 

0.114 

/ 

/ 

/ 

/ 

Rfcw_cond 

kW/(m2  K) 

1.561 

1.376 

1.373 

1.588 

1.658 

1.513 

1.64 

1.969 

Rfcw_supc 

kW/(m2  K) 

0.810 

0.858 

0.8892 

0.753 

0.8456 

0.9104 

0.9551 

0.5115 

Ac 

m2 

8.164 

7.975 

8.493 

8.531 

8.989 

8.776 

9.93 

11.57 

Ae 

m2 

10.24 

12.62 

13.24 

13.77 

13.43 

13.71 

14.21 

14.26 

Ank 

Dollar 

52,910 

60,821 

66,985 

70,479 

73,059 

75,299 

78,046 

81,092 

NPV 

Dollar 

-84,032 

24,582 

147,710 

309,068 

66,430 

64,857 

98,865 

88,035 

Appendix  B 

Evaluation  results  of  Case  Study  B. 

Single-loop  sub-ORC 

Single-loop 

sup-ORC 

Double-loop  ORC 

Units 

Pmax  2  MPa 

3  MPa 

6  MPa 

8  MPa 

Units  2  MPa 

4  MPa 

Units 

2  MPa 

4  MPa 

rhf  kg/s 

0.186  0.188 

0.1891 

0.1903 

rh/i 

kg/s 

0.03435  0.032 

ratiolcout 

% 

8.816 

8.884 

hr  kj/kg 

250.000  251.000 

253.8 

255.7 

mf  2 

kg/s 

1.185  1.181 

ratioIc2 

% 

15.080 

15.282 

h2  kj/kg 

306.700  306.600 

306.8 

307 

Pexpl 

kW 

21.450  22.270 

ratiolgfi 

% 

26.938 

21.868 
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Appendix  B  (continued) 


Single-loop  sub-ORC 

Single-loop 

sup-ORC 

Double-loop  ORC 

Units 

Pmax 

2  MPa 

3  MPa 

6  MPa 

8  MPa 

Units 

2  MPa 

4  MPa 

Units 

2  MPa 

4  MPa 

h3 

kJ/kg 

344.900  345.000  / 

/ 

kW  15.230 

15.170 

ratiolgp 

% 

4.758 

8.164 

kj/kg 

485.100  488.800  / 

/ 

m 

%  20.110 

22.370 

ratioIpumpl 

% 

0.031 

0.060 

h4 

kJ/kg 

973.100 

970.700 

963.7 

959.3 

m 

%  8.270 

7.939 

ratiolpumpi 

% 

0.163 

0.165 

h5 

kj/kg 

905.400 

891.400 

865.6 

853.8 

Qgi 

kW  106.700 

99.540 

ratioIeXpi 

% 

15.640 

17.416 

h6 

kj/kg 

431.500 

431.500  431.5 

431.5 

kW  17.040 

24.190 

ratioIex p2 

X 

12.418 

12.559 

h7 

kj/kg 

248.300 

248.300 

248.3 

248.3 

Qfw 

kW  167.100 

166.900 

ratioljp 

% 

11.274 

11.426 

Pexp 

kW 

12.610 

14.860 

18.56 

20.08 

hi 

kj/kg  315.8 

318.4 

ratioIfv 2 

% 

3.861 

2.965 

PfP 

kW 

0.316 

0.497 

1.04 

1.407 

h2 

kj/kg  908.7 

1087 

HL-Kgfievp 

l<W/(m1 2 3 4 5 6 7 8 9 10 11 12 

;K)  0.1137 

0.1148 

Qg 

kW 

124.200 

124.500 

124.2 

124.2 

hi  g 

kj/kg  2799 

2801 

HWWe 

kW/(m2 

!K)  0.09875 

0.09979 

Qjw 

kW 

10.560 

10.440 

10.02 

9.769 

h3 

kj/kg  3421 

3397 

HL-Kgpsupn 

kW/(m2 

!K)  0.093 

0.094 

/system 

kW 

48.180 

46.150 

43.19 

41.76 

h4 

kj/kg  2794 

2703 

LL-Ksupc  2 

l<W/(m2 

!  K)  1.482 

1.479 

ratiolcond 

% 

73.300 

73.790 

73.32 

73.42 

h3 

kj/kg  313.3 

313.3 

LL-Kcond2 

l<W/(m2 

;K)  2.736 

2.731 

ratioljp 

% 

0.126 

0.206 

0.522 

0.679 

h5s 

kj/kg  2634 

2634 

LL-Kffj, 

kW/(m2 

;K)  0.1543 

0.1463 

ratiolgf 

% 

18.330 

16.490 

13.7 

12.80 

h6 

kj/kg  248.6 

248.6 

LL-Kjjl 

kW/(m2 

!K)  0.1694 

0.161 

ratioljwf 

% 

1.379 

1.409 

2.07 

1.57 

h7 

kj/kg  389.7 

390 

LL-%_evp 

kW/(m2 

!K)  0.09163 

0.0939 

ratiolt 

% 

4.827 

6.039 

8.32 

9.44 

h7s 

kj/kg  456.9 

453.6 

LL-Kgfsuph 

l<W/(m2 

!  K)  0.09692 

0.09852 

rationlcoh 

a  % 

2.040 

2.065 

2.09 

2.092 

h8 

kj/kg  461.6 

455.5 

LL-KjwfH 

l<W/(m2 

;K)  1.909 

1.909 

Kfcw_cond 

kW/(m2 

!K)  0.840 

0.843 

1.042 

1.046 

h9 

kj/kg  291.9 

291.9 

LL-Kjwfj. 

l<W/(m2 

;K)  1.359 

1.356 

Kfcw  supc 

kW/(m2 

:  K)  0.734 

0.731 

0.7304 

0.7284 

hw 

kj/kg  453.7 

453.7 

HL_Ae 

m2 

11.010 

10.690 

Kgf  evp 

kW/(m2 

:  K)  0.101 

0.103 

/ 

/ 

hit 

kj/kg  476 

476 

LL_AC 

m2 

13.890 

13.810 

Kgf-pre 

kW/(m2 

:  K)  0.085 

0.086 

/ 

/ 

ftl2 

kj/kg  462.8 

462.8 

LL-Am 

m2 

71.430 

70.370 

Kgf_suph 

kW/(m2 

1  K)  0.103 

0.103 

/ 

/ 

hi  3 

kj/kg  248.3 

248.3 

LL-Agf 

m2 

5.072 

3.347 

Kgf 

kW/(m2 

!  K)  / 

/ 

0.07338 

0.07398 

hi  3s 

kj/kg  431.5 

431.5 

LL_Ajwf 

m2 

9.240 

9.234 

Kjwf 

kW/(m2 

:  K)  0.4963 

0.4974 

0.4983 

0.4991 

/system 

kW  84.16 

84.19 

Ank 

Dollar 

162,800 

161,800 

Ac 

m2 

4.765 

4.773 

4.09 

4.097 

NPV 

Dollar 

260,030 

275,023 

Ae 

m2 

22.350 

26.030 

48.1 

47.46 

Ajwf 

m2 

1.176 

1.169 

1.146 

1.132 

Ank 

Dollar 

75,016 

76,832 

86,082 

86,568 

NPV 

Dollar 

248,013 

299,966 

-44,953 

-34,006 
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